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Q^, Abstract. We study the asymptotics of the average number of squares 

■^ \ (or quadratic residues) in Z„ and Z* . Similar analyses are performed for cubes, 

pg ' square roots of and 1, and cube roots of and 1. 

(N' 

Let Tin denote the ring of integers modulo n, and let Z* denote the group (under 

f-H ' multiplication) of integers relatively prime to n. The number of elements in Z* is 

| 7 | ■ '^{f^)-i where y) is Euler's totient function. What is the average number of elements in 

r^ . Z*, given an arbitrary n? One way to answer this question is to apply the Selberg- 

c^ ! Delange method [U 12 IS| to the Dirichlet series 
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where G{s) is bounded in a half plane Re(s) > c for some c < 1. In fact, G{s) 
l/Cl-s + 1) in this case, and hence 

k " r(i) c(2) 

as A^ ^ oo. It follows by partial summation that 

^^^ ' 2C 2 7r2 

n<N ^^ ' 
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A more elementary proof of this asymptotic formula appears in 4]. Since the Selberg- 
Delange method will be used throughout this paper, we choose to illustrate its appli- 
cation in this simple setting. 

Many other questions can be asked for arbitrary n: 

• What is the average number of solutions of x^ = 1 in Z* ? 

• What is the average number of solutions of x^ = 1 in Z* ? 

• What is the average number of solutions of x^ = in Z„? 

• What is the average number of solutions of x^ = in Z„? 

• What is the average number of images of the map y ^^ y'^ m either Z,„ or Z* ? 

• What is the average number of images of the map y t— > y^ in either Z„ or Z* ? 

Although the answers require only straightforward use of standard techniques, 
they do not seem to be explicitly given in the literature. We make no claim of 
originality: Our purpose is only to collect results in one place and to document relevant 
numerical techniques. 

1. Number Theory 
1.1. Selberg-Delange Method. Let F{s) = ^^i a(^)^~^ be a Dirichlet series 
with positive coefficients and with the property that G{s) = F{s) ■ ({s)~^ can be 
analytically continued and is bounded over Re(s) > c, for some c < 1 and some 
z e C. Then 

n<N ^ ' 

as N ^ oo. More terms of the asymptotic expansion are possible, as is an accurate 
estimate of the error, but we omit these details for brevity's sake. 

A generalization of this method is required for our work involving averages over 
arithmetic progressions. Let x denote the principal character modulo k = q"^, where 
g is a prime and m > 1. Here we examine 

^ X{n)a{n) 
^xis) = 2^ = G^{s) ■ L^{s) 

n=l 

where 

n=l ^ ^ 
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is the L-series corresponding to x- Assuming a{n) is a multiplicative function and 
g f £, it follows that 

n<N, ^ ' 

n=£ mod k 

as N ^ CO and, further, that 

r=l 



oAn^^M.^^-^] «a,. 



In our examples, q will be either 2 or 3 and the bracketed infinite series will always 
collapse to a closed-form expression. Rather than directly employing the formula for 
G^{1), however, we prefer instead to deduce F^{s) (and hence Gy^{s)) from F{s) on 
basic principles. 

1.2. Square Roots of Unity. The number a{n) of solutions of x^ = 1 in Z* is 



2^(n)-i ifn = 2,6mod8, 

2'^W if n = 1,3, 4, 5, 7 mod 8, 

2<^(n)+i ifn = 0mod8 



where uj{n) denotes the number of distinct prime factors of n. It is well-known that 



and hence 



°° 2^W / oo ^ 

n=l "' p \ r=l ^ 



y- 2'^W ^ -L^N ■\nN = ^N ■ InN. 

n<N ^^ ' 



2 



We need to generalize this asymptotic formula to arithmetic progressions n = £mod /c, 
where fc = 2*" for m > 1 and 21^. It can be shown that 



2^(") ^ T\( 2 \ 



n=l.nioAk p>2 



fc V 2^-iy C(2s) ^^ ^ ^ ^ 



Squares and Cubes Modulo n 
as s — s> 1, and thus 

^ r 2 A;7r2 

n<N, ^ ' 

n=£niod A: 

The cases (/c, €) = (8, 1), (8, 3), (8, 5) and (8, 7) follow immediately. The case (A;, 
(8,4) proceeds from the case {k,i) = (2, 1): 

^ "S^ 2'^*^"^ = ^ \^ 2'^^")+^ > - 2 — = — 

N-lnN ^ AT -In AT ^ 4 vr^ 7r2' 

n<N, inKN, 

n=4mod8 ?i=lmod2 

The case {k,i) = (8,2) proceeds from the case {k,i) = (4, 1): 

^ >p 2'^{n) ^ 1 V^ 2'^(")+i > i 2 — = — 

A^-lnA^ ^ A^-lnA^ ^ 2 tt^ tt^ 

n<Ar, 2n<Af, 

71=2 mod 8 n=lmod4 



and (8, 6) likewise proceeds from (4, 1). By everything proved thus far, we have 






n<Af, 
n=0mod8 



Therefore 



5^ „(„).(!. 2. i-H. (4.^ + 1) +2. i)iV. in iV^AiV.lniV. 

n<Af ^ ^ ^ ^ 

It is interesting that (8, 2) and (8, 6) balance perfectly against (8, 0) so that the mean 
value of a{n) is asymptotically equivalent to the mean value of 2^^'^\ 

1.3. Cube Roots of Unity. The number a{n) of solutions of x^ = 1 in Z* is j7] 

J 3^{n) if n = 1, 2, 3, 4, 5, 6, 7, 8 mod9, 

aw- I 3^(n)+i ifn = 0mod9 

where ijj{n) denotes the number of distinct primes of the form 3/;; + 1 dividing n: 

ifp = 3or]9 = 2 mod 3, 



^1 if p = lmod3. 
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First, note that 



E— - n i+E- ■ n iH-3E^ 

Z_^ y^s J. J. \ Z-^ ipTS j J. J. \ A^ prs 

n=l p=3 or \ r=l / p=lmod3 \ r=l 

p=2mod3 



n (i+^)- n 



1 ' 



P'^ — 1/ -'--'■ v p'^ — 1 

p=3 or ^ ■' ' p=lmod3 

p=2mod3 

= n i^-^Y- n 64)"7-i)(^ 



p=3or ^ ^ ^ p=lmod3 ^ "^^ ^ /'/V /^V/^'y 

p=2mod3 

and IE] 

Si n 64)"'^-(-i)-# n 64^"' 



. p^ J 27r -'■-'- V n^ 

p=lmod3 p=lmod3 



hence 

V 3'^(") ~ ^i^N ■ InN = C-N-lnN 

n<N 

where 



r(2) 



C = — n fl-- -)= —(0.9410349413195354517900322...). 

p=lmod3 ^ ^^^ '^ 

We need to generahze this asymptotic formula to arithmetic progressions n = £ mod fc, 
where fc = 3™ for m > 1 and 3 f ^. It can be shown that 

3'^(") 1 T-r /. 1 \ T-r / 3 



3^ ^ -j-j- / 1 \ -j-r 



(j9(A;) -'--'- V p'^ — 1 J -'--'- \ p'^ — 1 

ji=^modfc p=2mod3 p=lmod3 

3 /. 1 ^ -^ 



7.^ 1 + 



as s — > 1, and thus 



y- 3-(n) _ ^illliv. In iV=^iV- In AT. 

^^ r 2 k 

n<N, ^ ' 



n=^modfc 
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The cases (fc, £) = (9, 1), (9, 2), (9, 4), (9, 5), (9, 7) and (9, 8) follow immediately. The 
case (9, 3) proceeds from the case (3, 1): 



i \^ Q-^Cn) = ^ "^ 3' 

, In Af ^^ N .\^N ^ 



N -hiN ^^ " N-\aN ^^ ~' 3' 3 9 

n<N, 3n<N, 

n=3mod9 n=lmod3 

and (9, 6), (9, 0) likewise proceed from (3, 2), (3, 0). Therefore 

Y^ a{n) ~ ( 8 ■ - + 3 • - W ■ InA^ = —C -N -InN = {0.317. ..)N ■ InN. 

n<N ^ -^ 

Unlike earlier, the mean value of a{n) is asymptotically greater than the mean value of 
3'^("), Our estimate improves upon Cloitre [7j, who gave (0.4...)A^-ln(A^) on empirical 
grounds. 

1.4. Squares in Z*. Let a{n) be as defined in section |1.2j . The number of 
squares, that is, the cardinality of images under the map y i-^ y^ in Z* , is 



b[n) 



a[n) 



( ifi{n) if^^2,6mod8, 

^Vy if n = l,3,4,5,7modi 

2u;(n)+l 



u, >_,, i, ^, 



if n = 0mod8. 



First, note that 



n=l 



7lP 



hence 



y ^H ^ 4^ AT . (lniV)-i/2 = C-N- (lniV)-V2 
Z^ n2'^H r 1/2 ^ ' ^ ' 



J(l/2) 

n<N ^ ' ' 



where 



C = ^TTfl + — ifl--) = ^(0.8121057111631225117062509.. 
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It follows by partial summation that 

y 44 ~^.iv^. (In ivrv^. 

n<N 

We need to generalize this asymptotic formula to arithmetic progressions n = £mod k, 
where k = 2"^ for m > 1 and 2 \ i. It can be shown that 

(/9(n) 1 T-r / p — 1 






^s+i2'^(") m(A;) J-J- \ 2n(n* - 11 



K'^3(?^)'''<^'-^<-""^^'''<^'-^<^'"^ 



as s ^ 1, and thus 

^ n2'^H r 1/2 ^ ' 5k ^ ' 

n<N, ^ ' > 

n=l mod k 



or 






2<^(n) 5^ 

n<Af, 
n=£ mod fc 



The cases (/c, £) = (8, 1), (8, 3), (8, 5) and (8, 7) follow immediately. The case (/c, 
(8,4) proceeds from the case {k,(^) = (2, 1): 

{InNy/^ y. ^{n)_ _ {InNy/^ y. 2y^{n) ^ 1_ ±^ _ ^ 

Jp Z^ 2^ ~ 'pp Z^ 2^(")+i ^ 16 ' To ~ 40' 



n<N, 4n<N, 

n=4mod8 ?i=lmod2 



The case {k,i) = (8,2) proceeds from the case {k,i) = (4, 1): 

(lniV)i/2 v'(n) _ (lnA^)i/2 ^ ^{n) ^1 1 4 C 

Jp Z^ 2'^(") ~ JP ^ 2'^(")+i ' 4 ' 2 ' 20 ~ 40 



n<N, 2n<N, 

n=2mod8 n=lmod4 



and (8,6) likewise proceeds from (4, 1). By everything proved thus far, we have 

(lnAr)V2 ^(^ ,^_A £_£_2 £-£ 

Jp Z^ 2^in) 'Y 'To~40~ ■ 40 ~ 40' 



n<N, 
n=0 mod 8 



Therefore 

yb{n) ~ (2.2--+(4-- + -]+---]N'-{\nN)-'/' 
^^^ ' V 40 V 10 40; 2 40; ^ ^ 



43 

—C-N^- (In Ar)-V2 _ (0.246.. .)iV' ■ {\nN)-^'^. 

80 
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1.5. Cubes in Z*. Let a{n) be as defined in section jl.3j . The number of cubes, 
tliat is, tlie cardinality of images under the map y ^^ y^ m. Z* , is [TU] 






V{n) J ^iiH if r^ = l,2,3,4,5,6,7,8mod9. 



3(:)(n)+i 
First, note that 



a[n] I Vin) -r _n ^n 

if ra = 0mod9. 



p=3 or 
p=2mod3 



V ^W ^ TT i + y-P^ . TT 1 + iV^^ 

n=l p=3 or \ r=l " / p=lmod3 \ r=l " 



TT fi + -^^^ . TT fi + P-^ \ = G{s) ■ C(.)^/^ 



p=3 or ^ ^ '' p=lmod3 

p=2mod3 



hence 

y ^H ^ 4^ AT . (lniV)-V3 = C.N. (lniV)-V3 



nS-^W r(2/3)' 



where 



p=2mod3 

^ -(0.9477556177621765519078142...). 



r(2/3) 

It follows by partial summation that 

y^~^.iV^.(lniV)-V3. 

Z^ 3<^(n) 2 

n<Ar 

We need to generalize this asymptotic formula to arithmetic progressions n = i mod k, 
where fc = 3*" for m > 1 and 3 \ i. It can be shown that 

filT-) 1 TT f -. P " ^ \ TT A p — 1 



V ^^^^^ — T[ (1+ ^"^ V TT (1 + 



._ i3'^(") (z)(A;) J- J- V P(P'-I)/ ■'■■'■ V 359(n*-r 

n^fmodfc ^^ ^ p=2mod3 ^ ^^^ '' p=lmod3 ^ ^^^ ' 
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as s — s> 1, and thus 






nS-^W r(2/3) ' ' 8k 

n<N, ^ ' ' 

n=l mod k 



or 






n<N, 
n=£Taodk 

The cases {k, i) = (9, 1), (9, 2), (9, 4), (9, 5), (9, 7) and (9, 8) follow immediately. The 
case (9, 3) proceeds from the case (3, 1): 

(lniV)^/3 ^ !f{n) _ (lnAr)i/3 2^(n) ^1 _9_^ _ £ 

Jp Z^ S'^H ~ Jp 2^ 3'^(") '9' '48 ~ 24 

n<N, 3n<N, 

n=3mod9 n=lniod3 

and (9,6) likewise proceeds from (3,2). By everything proved thus far, we have 

(lniV)V3 ^^^_Q £_2 - = - 

N^ 2^ 3<^(n) '2 ' 16 '24 24' 

n<N, 
n=0mod9 

Therefore 

n<N ^ ^ ' ' 

= YLc-N''- (lnAr)-V3 _ (0.330.. .)Ar2 ■ (lnAr)-V3. 
36 

1.6. Square Roots of Nullity. The number a{n) of solutions of x^ = in Z„ is 
a multiplicative function of n, with a(p'") = pL''/2J^ ^];ius [TTI[T2] 

Eo(n) T-r /_ 1 p J) r? r? rr' p^ 

-^^ = IT 1 + — + -^ + -^ + '^ + '^ + 'ir + '^ + - 
lY^ ^ .A. .X. V 7^ lyi^S ly-iOS ^4tS ^DiS iT-jOS /y-t I S 

n/ ]^ p ]^ „ „2 ]^ 2 3 ]^ „3 \ 

1 + — 1 + -^ + ^ + ^ + —^ + ^ + —^ + ■ ■ ■ ] 
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10 



and lims-+i C(2s — 1) ■ (s — 1) = 1/2, hence 

G(l) 



n<N 



a(n) 



r(2) 



N -InN = —N-\nN. 



ir^ 



1.7. Cube Roots of Nullity. The number a{n) of solutions of x^ = in Z„ is a 
muhiphcative function of n, with a(p^) = pL2f/3j^ thus ^H] 



E 

n=l 



a(n) 



77," 



+ 



n/ ]^ p p2 2 3 „4 4 5 pi 

IH h-^ + -^- \- — \- — \- — \- — \- — h — 
\ 7^^ T) T) 7^ T) Tl f) f) TY' 



,9s ' plOs 



n 






Y r^is 



pi 



1 



p- 



p- 



p 



,3s 



p 



,3s 



n(i-^ 



-1 



1 1 

1+-+ 



pS p. 



,2s-l 



C(3.-2)n(l + ^ + ^)=G(.)-C(^ 



We have 



lim TT ( 1 + — + ^-— 



's - 11 



i'-^c(^)-2C(2s-i)n(i+ps+p2s-ij 

mo 

p ^ 






, 1 - - 
p — I J \ p 



m(-^)0 



p(p + 1) 



2C(2) n (^ 



p(p + 1) 



and linis^i ({3s — 2) ■ (s — 1) = 1/3, hence 

G(l 



E «(^) ~ ^^ ■ (1^^)' = ^ ■ ^ ■ (1^^)' 



n<N 



r(3) 
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where jHl II3I 



p{p+l)J 12 



— (0.2867474284344787341078927...). 



1.8. Squares in Z„. The number h{n) of images under the map y ^^ y^ in Z„ is 
a multiplicative function of n, with [TU [121 CHI 



6(/ 



- (2"^-^ + 4) 

? 

3 (2^-^ + 5) 



if p = 2 and r = mod 2, 
if p = 2 and r = 1 mod 2, 
if p > 2 and r = mod 2, 



^ r ip'^^ + 2)0+1) if n > 2 and r = 1 mod 2 

2 n + 1 



and 



I 2(p + 



6(n 



6(2^ 



6(/ 



n=l \ r=l / p>2 \ r=l ^ 



The left-hand factor in F{s) simplifies to 

1 ^ 2^^-i + 4 1 ^ 2(^-^'-i)-i + 5 



i=l 



1 



4 



s+l 



+ 2^^ 



2 ■ 4^+1 - 7 



2 \^(4«+i- 1)(4^- 1) (4^+1-1) (4^-1) 

22S+1 _ 2^+1 _ 1 ^^ ^^ 2^+^ + 2 

2s+2(22«+i - 2^-1 - 1) 



1- 



2(2«+i + l)(2^+i-l) 



1 

2^ 



-1 



and the p^^ right-hand factor simplifies to 



1 + 



1 ^p^'+^+p + 2 



ttE 



^ 00 



p 



,(2i-i)+i 



2j9+l 



1 + 
1 



2(p+i)Z^ p(2i){^+i) ^2(p+l)4^ p(2i-i)(.+i) 
1 /p2s+3 ^ p2s+i ^ 2p2^ - 2p - 2 ^ ,+iP^'+^ + 2p2^+^ + p2« - 2p - 2 



2(P+1)V (p2«+2_l)(p2«_l) 

(p^+i + 2)(p-l) \ / _ l_- ' 

2(j9«+i + i)(j9«+i-i)y V r 



+ p' 



(p2^+2-l)(p2s_l) 
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We have 



F{s) = as) 1 + 



)2s+l os+1 



25+2(22'^+! - 2^-1 - 1) 



n^ 



(p^+i + 2)(p-l) 
2(p«+i + l)(p^+i-i; 



G(s)-C(s 



,1/2 



and hence 



where 



y M ^ G(l) _ ^1^^-1/2 = c-N-{\n iV)-V2 
^^ n r 1/2 ^ ^ ^ ^ 



n<N 



c 



17 


1 


16 


v^ 


17 


1 



n^ 



p2 + 2 



2(p2 + i)(p+i) 



1-i 
P 



-1/2 



16 v^ 

It follows by partial summation that 

C 



1.2569136102101885959492115...). 



J2 b{n) ~ - ■ Ar2 . (lnAr)-i/2 = (0.376.. .)iV2 ■ {\nNy'/\ 



n<N 



1.9. Cubes in Z„. The number b{n) of images under the map y \-^ y^ in Z„ is a 
multiplicative function of n, with 



r 1 

1 

13 



(3''+i + 10) 

{y+^ + 30) 
(3^+1 + 12) 



p2 + 



p + 1 



6(p'-) = <^ 



pz 



f 



p^ + p + 1 



(p''+2+p+l) 
^pr+2 _^p2 _^p-^ 



if p = 3 and r = mod 3, 

if p = 3 and r = 1 mod 3, 

if p = 3 and r = 2 mod 3, 

if J9 = 2 mod 3 and r = mod 3, 

if p = 2 mod 3 and r = 1 mod 3, 

if p = 2 mod 3 and r = 2 mod 3, 



1 



3(p2 +P+ 1^ 
, 3{p'^+p+l] 



{p 



r+2 



2p2 -|- 3p -|- 3) if p = 1 mod 3 and r = mod 3, 



(jf+'^ _j_ 2p2 _|_ 3p _|_ 2) if p = 1 mod 3 and r = 1 mod 3, 
(jf+'i _|_ 3p2 _|_ 2p -|_ 3) if p = 1 mod 3 and r = 2 mod 3 
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and 

n=l \ r=l / p=2 \ r=l ^ / p=l \ r=l ^ 

mod 3 mod 3 

The expressions for b{p^) follow from a conjecture by Wilson ^\] a proof for the case 
p = 2, r = mods was given by Wilmer & Schirokauer ^Hj- The left-hand factor in 
F{s) simplifies to 

1 /y. 33^+^ + 10 ^ 3(3i-2)+l ^ 30 ^ 3(3fc-l)+l ^ 1 2 
^ + 13 [1^ 3{3i){s+l) + 1^ 3(3i-2)(s+l) + 1^ 3(3fc-l)(s+l) 

\j=i i=i fe=i 

^ I . 7-2r-l ^ ^,,^, 9-2r-7 ,^,+1 3-2r-l 



(27^+1-1) (27^-1) (27''+i-l)(27^- 1) (27^+i - 1)(27^ - T 



. 2 (3-+^ + 1) \ ( l_ 

(3^+1 + 3(^+i)/2 + l)(3«+i - 3(^+i)/2 + l)(3«+i - 1); V 3^' 



-1 



The p* right-hand factor simplifies to 

'^p2_^p^llz^ p(3i)(3+l) +2^ p(3i-2)(s+l) +2^ p(3fc-l)(s+l) 

\ t — -L J — J- rC — _L 

1 / ^oS-j-O I f^oS~\-l. I ^OS j-,Z ^ 1 



+ p 
+ p 



2s+2 



j-.OS-i'O _i_ ^ds + z _i_ ^ds-f-i ^z ^ 1 



(p3«+3_l)(p3s_l 

s+1 



(j-jdS-|-4t I ^dS-|-Z I jyjS j^Z ly. 1 



(p3«+3_l)(p3s_lj 



^s+1 j_ ^\i^ _ ^\ \ / i \ ~i 



. (p-+^ + l)(p-l) \ f l_ 

(ps+l+ p{s+l)/2 + I ) (pS+1 _ p{s+l)/2 + 1 ) (pS+1 - 1 ) yl \^ p^ 



when p = 2 mod 3 and 



(oo oo oo \ 

V^ p3'+2+2p2+3p+3 I V^ p(3i-2)+2^3p2^3p^2 , V^ p(3fc-l) + 2^3p2^2p+3 ) 
/ . p(3i){s+l) -r 2_^ p(3i-2)(s+l) r 2_^ p(3fc-l)(s + l) I 

i=l i=l fc=l / 

1 /p3.+5 ^ 2p3«+2 + 3p3s+l _^ 3p3s _ 3p2 _ 3p _ 3 



3(p2+p+l)\^ (p3s+3_l)(p3s_l) 

2s+2P^'"^^ + 3p^'+^ + 3p3"+i + 2p3« - 3p2 - 3j9 - 3 

^ (p3.+3_l)(p3. _1) + 
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/ (2p^-+^ + 3p-+^ + 3)(p-l) Vl-lV 

when p = 1 mod 3. We have 

2 (3^+^ + 1) 

(3S+1 + 3(s+i)/2 + i)(3s+i _ 3(s+i)/2 + i)(3«+i _ 1) 

(p^+i + l)(p-l) 



•n6 



_ ^ (ps+l + p{s+l)/2 + 1) (pS+1 _ p(s+l)/2 + 1) (pS+1 _ 1) 

P — 1 

mods 
mods 

and hence 

E ^ - tWM^ ■ ^1" ^)"'' = C^ ■ iV ■ (In iV)-V3 



where 



n r(2/3)' 

n<N ^ ' ' 



12 1 /. i\-V3 



(^ — \ 

13r(2/3) V 3 

mods 

TT A 2p^ + 3p^ + 3 \ / _ n -^/^ 

11 V 3(p2+P+l)(p2-p+l)(p+l)y' V V 

mods 

1.4225831466986636811460982...). 



13r(2/3) 
It follows by partial summation that 



^ hin) ~ - . Ar2 . (lnAr)-V3 _ (0.484.. .)Ar2 ■ (lnAr)-V3. 



2 

n<Af 



We emphasize that this result is only conjectural. 
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1.10. Other Problems. The power of the Selberg-Delange method is evident 
(many deeper applications occur elsewhere in the literature). We merely mention 
that the number a{n) of solutions of x^ = — 1 in Z* satisfies 

^ ' 27r ' 

n<N 

in particular, x^ = — 1 has asymptotically far fewer solutions than x"^ = 1. Such 
asymmetry does not occur for x^ = ±1 (just replace x by —x). See other modu- 
lar polynomial equations at ^Hl and the enumeration of weakly primitive Dirichlet 
characters at [201 |2I1 • 

A more difficult exercise concerns the number b{n) of elements of Z„ that are 
both squares and cubes. li w = z^, then clearly w = [z^Y = [z'^Y- Conversely, if 
w = v^ = f^, then {uv~^Y — {u'^Y{v'^)^'^ = w^w~'^ = w. Hence b{n) is the same as 
the number of sixth-powers in Z„. Wilson's conjecture again provides expressions for 
b{p^), which in turn give formulas for F{s) and G{s). The details of this and other 
higher-power problems are left to someone else [22] ■ 

2. Numerical Techniques 
2.1. Prime Products. Here is a method for evaluating constants of the form 

C= U f{p) 

p=i mod k 

to high precision, where the product is taken over all primes of the form p = mk + i. 
Suppose that the function In / has asymptotic expansion 

as p -^ oo, where (c„, s„) are real numbers and 1 < S2 < ... < s„ < .... (Often Sn = n 
occurs.) Define the {k,iy^ prime zeta function 



PkAs) 



p=£ mod k 



for Re(s) > 1; it follows that 



lnC= J^ ln/(p) = J^c„Pfc,<?(s„). 

p=^modA; n>2 

Let pk/ denote the smallest prime of the form mk + £; clearly Pk,e{n) ~ ^/Pke ^s 
n ^ oo. Consequently, if the coefficients c„ are uniformly bounded, the convergence 
of the sum is fast (geometric). It hence remains to accurately compute the values 

PkA^n)- 
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2.2. Prime Zeta Functions. Let Re(s) > 1. The classical prime zeta function 
P(s) = Pifl{s) can be related to the classical zeta function by Euler's famous product: 






Pins) 



n 



Applying the Mobius inversion formula, we obtain [221 124 

/i(n) 






n 

n=l 



Since ln(^(ns) ~ 2~"'* as n ^ oo, only a few terms in this series are required to 
compute an accurate value of P{s). Also P{s) ~ — ln(s — 1) as s ^ 1^. These facts 
are useful in computing constants of the form Y\pf{p)- 

For constants of the form np=£mod3 /(^)' "^^ need -P3,i(s) and P3,2(s). To achieve 
this, it is necessary to introduce the two characters modulo 3: 

1 ifn=lmod3, ( 1 ifn = lmod3, 

Xo{^) = { 1 ifn = 2mod3, Xi{^) = { ~1 ifn = 2mod3, 

ifn = 0mod3, [ ifn = 0mod3, 

and their associated Dirichlet L-series: 

n=l 'I' <-> 

where ( (s, a) is the Hurwitz zeta-f unction. By well-known acceleration procedures, 
series of this nature can be evaluated to many decimal places. 
From the Euler product expressions 



p ^ 

^i(^) = n (i 



XoiP)\ 
P' J 


"- n 0--) 

ii \ psi 

p=lmod3 ^ ^ ^ 


1 

p=2mod3 ^ -^ ^ 


-1 


Xi{p)\ 

p' J 


"- n fi--) 

p=lmod3 ^ -^ ^ 


p=2mod3 ^ -^ ^ 


-1 



we obtain 



2^\Li(s)J 2 ^ ^ll-P^V ^ 2n+l 

^ J^V ''/ p=2mod3 ^ -^ ^ n=0 
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l^^/ Lo(g)Li(s) \ _ 1 >^ ^^/ l+p-^ \ _^ P3.i((2n + l)g) 
2"^^ ^o(2s) / 2 ^ "^U-P^V ^ 2n + l 

^ UV 7 / p=lmod3 ^ -^ ^ n=0 



and, again, by Mobius inversion. 






2^ 2n + l V^i((2^ + l)s 



n=0 



l>^/i(2r2 + l) f Lo{{2n + l)s)L,{{2n + l)s) \ 
'^^^^^'^-2^^^^Tr^H ^o((4n + 2).) J- 

Of course, Lo(s) = C(s)(l - 1/3") and Li(s) = 1 - 1/2" + 1/4" - 1/5" + ■ ■ ■ . Also, 

^3,1(5) i ln(s - 1) and ^3,2(3) 1 ln(s - 1) as s ^ 1+. 

Similar techniques involving characters modulo k can be used to compute con- 
stants of the form rip-^modfc /(p)? but for brevity's sake we do not discuss these here. 

2.3. A Simple Example. Let us compute the constant 

c= n (1 ^ )- n ( (p-^)(p+^) ] 

piidsV pip+^)J JLA ip + ^)p J 

that appears in section |1.3j . It is easy to establish that 

">^- E (ln(^)..n(l.H))^E^^..(n) 

p^lmodS V V-f^ ' / V -f^/ / n>2 

where c„ = 2" — 2 when n is odd and c„ = —2^ when ra is even. Since c„ = 0(2") and 
-f3,i(^) = 0(7^"), it is more efficient to compute directly the product up to a certain 
cutoff pc- For example, if we take pc = 31, we find 

3247695 ^ / 



3430336 J- J- V p(p + l) 

p=lmod3, ^ ^^^ ' 

p>31 



and consequently 



3247696^5,/ l_^_^_n 

3430336 ^ n\ '^ ' 7" 13" 19" 31" / 

n=2 ^ ^ 
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enjoys much faster convergence because P^^i{n) — ^ — -^ — -^ — -^ = 0(37 "). 
The first few terms of this series produce 



n 



2 
3 

4 
5 

10 
15 



C 



0.94(09438379523896292195206.. 
0.94103(87732177050567463275.. 
0.941034(8096648041499806620.. 
0.94103494(70255355752383278.. 
0.94103494131953(43277214763.. 
0.941034941319535451790(3566.. 



and only 15 terms are necessary to obtain 20 correct decimal places. 
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